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Thermalization of isolated many-body systems is demonstrated by generalizing an approach orig¬ 
inally due to von Neumann: For arbitrary initial states with a macroscopically well-defined energy, 
quantum mechanical expectation values become indistinguishable from the corresponding micro- 
canonical expectation values for the overwhelming majority of all sufficiently late times. As in von 
Neumann’s work, the eigenvectors of the Hamiltonian and of the considered observable are required 
to not exhibit any specially tailored (untypical) orientation relative to each other. But all of von 
Neumann’s further assumptions about the admitted observables are abandoned. 
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The universal and irreversible tendency of nonequilib¬ 
rium states towards thermal equilibrium is an everyday 
experience in the macroscopic world, but in spite of more 
than a century of theoretical efforts, it has still not been 
satisfactorily reconciled with the basic laws of physics, 
which govern the microscopic world, and which are fun¬ 
damentally reversible [l|. The first quantum mechanical 
exploration of this problem is due to von Neumann Q, 
was unfortunately misunderstood for decades, but has 
recently been rehabilitated in a very enlightening com¬ 
mentary by Goldstein, Lebowitz, Tumulka, and Zanghi 
Q. A major remaining bottleneck of von Neumann’s ap¬ 
proach is his notion of “macro-observer” or “macroscopic 
measurement”, stipulating that all relevant observables 
can be approximated by commuting Hermitian opera¬ 
tors with very high-dimensional common eigenspaces Q . 
As an alternative, Goldstein et al. i3 suggested to 
consider “macroscopic observables” with the additional 
property (excluded in von Neumann’s original treatment) 
that one of those eigenspaces is overwhelmingly large 
compared to all the others. In our present work, all such 
restrictions with respect to the considered observables are 
abandoned. 

As in Refs. M, we consider an isolated many-body 
system, whose energy E is known up to an uncertainty 
5E, which is small on the macroscopic but large on the 
microscopic scale. The system is modeled by a Hamilto¬ 
nian El with eigenvalues En and eigenvectors |n), n £ N. 
System states (pure or mixed) are described by den¬ 
sity operators p, evolving in time according to the usual 
Liouville-vonNeumann equation p[t) = i[p(t), E[]/h. Ob¬ 
servables are modeled by Hermitian operators A with 
expectation values {A)p := Tr{pA}. The preset energy 
interval [E, E -\- 5E] defines an energy shell, namely the 
Hilbert space EL spanned by all \n) with En £ [E, E-\-5E]. 
Without loss of generality, we assume that the corre¬ 
sponding labels are n = 1,2,... ,D. For a macroscopic 
system with, say, / ~ 10^^ degrees of freedom^he di¬ 
mensionality I? of "H is exponentially large in / [3| , sym¬ 
bolically indicated as 

( 1 ) 

By definition, the probability to encounter a system 
energy outside [E, E + 5E] is negligibly small, and is 


henceforth idealized as being strictly zero. As a con¬ 
sequence, the diagonal matrix elements (“level popula¬ 
tions”) pnn '■= {'n\p\n) vanish for all n > D, imply¬ 
ing with Cauchy-Schwarz’s inequality that pmn = 0 if 
m > D or n > D. Denoting by P the projector onto 
EL, the projection (or restriction) of A onto EL takes the 
form A := PAP and analogously El := PHP etc. (note 
that p = p). It readily follows that Tr{pA} = Tr{pA} 
and that H yields the same time evolution of p{t) as 
H. Hence we can and will restrict ourselves to the en¬ 
ergy shell EL from now on, but, for convenience, omit 
the tilde symbols. Accordingly, P becomes the iden¬ 
tity operator on EL and the microcanonical density op¬ 
erator follows as pmc ■= P/D with expectation values 
(^) me ■ — Tr{p 

me 

The problem of thermalization is to show that (A)pp) 
evolves towards (A)mc for arbitrary (possibly far from 
equilibrium) initial conditions /9(0) : EL ^ EL. It is well 
known that this is impossible without additional assump¬ 
tions on H and A. With respect to H, we adopt von 
Neumann’s assumption 3 that the energy differences 
Em — En are finite and mutually different for all pairs 
m ^ n. Excluding nongeneric cases with additional con¬ 
served quantities (besides H), e.g. due to (perfect) sym¬ 
metries or noninteracting subsystems, the validity of this 
assumption is by now commonly accepted [MU. More¬ 
over, one expects that even considerably weaker assump¬ 
tions will do [T^4T5l| . 

Denoting by Umax and amin the largest and smallest 
among the D eigenvalues of A, the range of A is defined 
as Aa '.= Umax — cimin- Furthermore, in any real (or 
numerical) experiment, {A) p can be determined only with 
some finite accuracy 5A. In practice, we thus can focus 
on measurements which yield at most, say, 20 relevant 
digits, i.e. 

SA > Aa 10“^° . (2) 

The eigenvectors of H and of A are related by some 
unitary basis transformation U. A key point of von Neu¬ 
mann’s approach is the assumption that these two eigen- 
bases do not exhibit any “special orientation” relative to 
each other i[a, i.e. , the actual U is “typical” Q among 
all possible unitary transformations U : EL ^ EL in the 
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following sense: If a certain property can be shown to 
hold for the vast majority of C/’s (uniformly distributed 
according to the Haar measure [2^ ll^ ) ^ then this prop¬ 
erty is supposed to hold for the actual U as well, unless 
there are special reasons to the contrary. Denoting by 
IJLu{X) the fraction (normalized measure) of all C/’s ex¬ 
hibiting a certain property X, a /ic/(X) value close to 
unity (zero) is thus assumed to generically imply (ex¬ 
clude) property X for the actual system. While a more 
rigorous justification is clearly very difficult, intuitively 
such a “typicality” argument is very convincing: If we 
imagine A as fixed and H as arising by randomly sam¬ 
pling its eigenvectors via U [1, [H, [T3] , the argument is 
essentially tantamount to the common lore of random 
matrix theory Q, which is well known to be extremely 
successful in practice [l 8 j |. In particular, iiu{X) may be 
formally viewed as the probability of observing property 
X for a randomly sampled H (or C/), however keeping 
in mind - exactly as in random matrix theory - that 
there is no random sampling procedure in the real physi¬ 
cal problem under consideration ii- In passing we note 
that von Neumann actually adopted the complementary 
viewpoint of considering H as fixed while varying the 
eigenvectors of A 0 . 

By exploiting the above mentioned common assump¬ 
tions about the ener gy e igenvalues [ 1 , 0 -[ll| , one can 
infer [Mil (see also [T^) that the quantity 

a2(/) := [(kl),(,) - (A),]" (3) 

satisfies the relations 

D 

= XI l/’mn(0)p < max|^m„p , (4) 

m^n 

m^n 

where Amn ■= (m-|^|n), Pmn(O) := (m|p( 0 )|n), and the 
overbar indicates an average over all times i > 0. In par¬ 
ticular, p := p[t) is an auxiliary density operator with 
matrix elements pmn = ^mnPnnifi), sometimes named di¬ 
agonal or generalized Gibbs ensemble [l^ . The so-called 
eigenstate thermalization hypothesis (ETH) conjectures 
that for a many-body system with f 1 degrees of free¬ 
dom, typical off-diagonal elements A^n in O are expo- 
nentially small in / d, ■ Within our present gen¬ 

eralization of von Neumann’s approach, we can actually 
prove that even their maximum on the right-hand side in 
dH) is typically so small that 

pu (<T^(/) > e) < 4 +2\nD^ (5) 

for any e > 0. Besides O, the key jMredient in deriv¬ 
ing this result is Levy’s lemma (see and further 

references therein), stating that 

Prob(|g(((i) - (g)| > e) <2 exp (0 

for randomly and uniformly distributed points (j) on the 
d-dimensional unit sphere S'’* C and any Lipschitz 


continuous function g : > K with Lipschitz constant 77 

and mean value (g). Furthermore, any normalized |0) S 
TL of the form can be represented (via the 

real and imaginary parts of the Cn’s) as a point 4> on the 

a ) — l)-dimensional unit sphere. Finally, one can show 
that g{(j)) := ((/)|A|(^) is Lipschitz continuous with 
r] = Aa and (g) = (A) me Observing that randomizing cj) 
is equivalent to randomizing U, we thus obtain 



The remaining task is to connect this result for ((/)|A|^) 
with the maximal |Am,ra| in (j4|). The details are rather 
straightforward but tedious and thus provided as sup¬ 
plemental material in [T^. As an aside, it follows that 
von Neumann’s main technical achievement (Appendix 
of i) , a s well as its further improvement by Pauli and 
Fierz [^As in fact quite closely related to Levy’s lemma 
(see also 

Equation represents the first main result of our 
Letter. By choosing, e.g., e = D~^/‘^6A^ in dS]), it fol¬ 
lows with dD and m that the time-averaged variance 
from dS]) remains extremely much smaller than SA^ for 
“almost all” U: The fraction of the exceptional U’s is 
an unimaginably small number of the order of 10 “''’ with 
X ~ / Ri 10^^. Furthermore, the mere existence 

of the infinite time average in dH) implies that a similar 
estimate must also apply to averages of cr^ {t) over finite 
time intervals [0, T] with sufficiently large T BlMliil. 
Finally, the smallness of the latter time average implies 
(obviously or by Markov’s inequality) that the averaged 
quantity itself must be exceedingly small for most 
times/G [0,r] @,0,113, in. For example, for our above 
choice e = D~^/^5A^ and sufficiently large T, all the “bad 
times” t G [0,T] with |(A)pp) — {A)p| > 5A add up to a 
set, whose Lebesgue measure is smaller by (at least) a 
factor of the order k, (cf. dB)) than the mea¬ 

sure of all / G [0,T]. Altogether, we thus can conclude 
that for the overwhelming majority of f/’s, the differ¬ 
ence (A)p(() — {A)p remains below the resolution limit 5A 
for the vast majority of times t contained in any suffi¬ 
ciently large time interval [0,T]. The same conclusion 
carries over to our actual Hamiltonian H and observ¬ 
able A, given their eigenbases are related by a “typical” 
transformation U as discussed above. To establish quan¬ 
titative bounds for T is a subject of considerable current 
interest [sl.[l3. [23 - l3ll| . but goes beyond our present scope. 

The salient point is that ® holds independently of the 
initial condition p(0). Once a pair H, A with a typical 
U is given, the above implications of dS]) thus apply to 
any p(0): No matter how far from equilibrium the system 
starts out, for almost all sufficiently late times it behaves 
practically as if it were in the state p. Such an apparent 
convergence towards a stea^ state has been denoted as 
equilibration, e.g., in Refs. |9l4l5l|. 

To demonstrate thermalization, we still have to show 
that the difference between {A)p and (A)mc is negligibly 
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small. Recalling the definitions of these two expectation 
values, one readily sees that 

D 

B := {A)p — (A)iiic = Prm(O) [Ann ~ (^)mc] (8) 

n—1 

and hence that 

\B\ < ma.x\Ann - (^)mc| ■ (9) 

n 

Similarly as above Eq. one part (actually the better 
known part) of ETH consists in the conjecture that typi¬ 
cal differences Ann — (^)mc are exponentially small in / 
[^, [20l - l^ . Within our present framework, we can prove 
that even their maximum in @ is typically so small that 

/i(7 (1^1 > e) < 2 expj-^^-|-lnD| (10) 

for any e > 0. This represents our second main result, 
whose derivation from © is quite obvious and is pro¬ 
vided in full detail as supplemental material in [l^ . Once 
again, it is crucial to note that (fTOll is independent of p 
(and thus of p(0)): Given a pair H, A with a typical U, it 
follows from ® and (|T0)l that the difference {A)p — {A)nic 
remains way below the resolution limit 5A for any p (or 
p(0)). Einally, upon considering A as fixed and H as 
a random matrix (see above) we can conclude that von 
Neumann’s approach H in fact anticipates the verifica¬ 
tion of ETH from Ref. [21| within a random matrix the¬ 
oretical framework, see also i[i3- 

One readily sees that the measure of all C/’s which give 
simultaneously rise to both equilibration, as discussed in 
the paragraph below 0, and negligibly small B values 
according to (fTUl) is still extremely close to unity. Hence 
thermalization follows for any given pair H, A with a 
generic relative orientation of the eigenbases, no mat¬ 
ter how the initial condition p(0) is chosen. Along the 
same lines, one can infer the simultaneous thermaliza¬ 
tion of several (not necessarily commuting) observables 
[iSlal, as long as their number remains “reasonable” 
(e.g., smaller than D). 

Similarly as for C/, let us now denote by V the unitary 
basis transformation between the eigenvectors of the den¬ 
sity operator p(0) and those of H. Likewise, pv{X) now 
represents the fraction (normalized measure) of all uni¬ 
tary transformations V : Tl ^ TL which exhibit a certain 
property X. Furthermore, the usual von Neumann en¬ 
tropy is defined as S'[p] := —fcsTrjplnp} and satisfies 
0 < 5'[p] < ^[pnic] = fcs InD. Hence the entropy range is 
As = ksAiD and, similarly as in ([2]), experimentally re¬ 
solvable entropy differences 5S can be assumed to satisfy 
5S > qAs = qkslnD for some small but still “reason¬ 
able” q value. It follows that p(0) entails a p with the 
properties that ^[pmc] ~ S'))?] > 0 and, as demonstrated 
in detail in the supplemental material [l^ . 

( 11 ) 


for any s > 0. This is our third main result. By choosing 
s = 6S and recalling that SS > qkBlnD (see above), it 
implies with o that the entropy of the diagonal ensem¬ 
ble p, towards which the “true” p{t) seems to equilibrate, 
differs from the microcanonical entropy only by an un- 
measurably small amount for a generic p(0), i.e. one 
without a specially tailored orientation of its eigenbasis 
relative to that of H. We remark that already von Neu¬ 
mann demonstrated a somewhat similar, so-called H the¬ 
orem [^, however, for a differently defined entropy, whose 
physical relevance has been questioned, e.g., in Ref. 0. 
Further related but different results about entropies of 
diagonal ensembles are also due to (3^ . 

As shown in Refs. [M, [H, an alternative up¬ 
per bound for the left-hand side of o is given by 
(Aj^/4) Tr{p^}. Eor the latter factor, Tr{p^}, a similar 
relation as in CB) is derived in the supplemental mate¬ 
rial 0 , yielding 

MV > e) < Al/(2ei^) (12) 

for any e > 0. By analogous arguments as in the discus¬ 
sion of (5) in the paragraph below ©, this amounts to 
an alternative demonstration of equilibration did. 
But in contrast to which applies to arbitrary p(0), 
provided the relative eigenbasis orientation of H and A 
is generic, the present findings now apply to arbitrary 
observables A, provided the eigenbases of H and p(0) are 
in a generic constellation. 

Finally, let us denote by W the unitary basis transfor¬ 
mation between the eigenvectors of the density operator 
p(0) and those of A, and consider 

:= [(Al(t) - (Ajmcf ■ (13) 

Similarly as before, we now can show (see supplemental 
material 0 ) that 

MW (o-mc(0 > e) < ^aK^D) (14) 

for arbitrary t and e > 0, and that 

it < Al/{eD) (15) 

for arbitrary ti < t 2 . Note that while (11^1) is a t inde¬ 
pendent upper bound for the measure of all VE’s with 
o'mc(^) ^ fdis does not imply that the set of those W’s 
is /-independent. An analogous caveat applies to (1151) . 
Yet another crucial point is that da and da are valid 
for completely arbitrary (even time-dependent) Hamilto¬ 
nians H -.n^n ElMllii. 

A first remarkable implication of (iMl) and (fTKll fol¬ 
lows by considermg A as “given” (arbitrary but fixed): 
Namely, “most” p(0) then yield practically the same 
expectation value {A)p(^t) as pmc for any arbitrary but 
fixed time point /, but also for practically all times t 
within an arbitrary but fixed time interval [ti , t 2 ] (see 
paragraph below ©)■ Put differently, nonequilibrium 


MV (S'[p me ] - -Sip] >s)< ks/s 
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expectation values are “untypical” (even for t = 0), they 
require very special orientations W of the eigenbasis of 
/o(0) relative to that of A. In particular, for pure states 
p(0) = !'(/’) ('01 we recover the quintessence of so-called 
canonical typicality and related phenomena [3^ . 113- 
(see also [Ibl. l4(1|l. 

Conyersely, when considering p(0) as “giyen”, it fol¬ 
lows from Cl and m that “most” [s^ measurement 
deyices A cannot distinguish p{t) from at any arbi¬ 
trary but fixed time point t, or for practically all times t 
within an arbitrary but fixed time interval [fi, ty]. This 
is the viewpoint adopted, e.g., in Refs. [H, |4l|, but now 
formulated within our present generalization of von Neu¬ 
mann’s original approach (see also 0). 

The fact that a nonequilibrium value of (^)p(o) requires 
an untypical pair p(0), A implies Q that conclusions re¬ 
garding thermalization can be drawn only from results 
concerning all orientations of p(0) relative to A, as it is 
the case in von Neumann’s approach (see below © and 
m), but not from results concerning most orientations, 
as in the above generalization ([T^ and p5l) of the ap¬ 
proach from Refs. [H, [H, [13, l4l| . In other words, it is 
not right to say that von Neumann’s approach is inad¬ 
equate to investigate thermalization since the same ap¬ 
plies to the approach from Refs. [H, m, 113, El . Rather, 
the two approaches are fundamentally different: One re¬ 
quires a generic eigenbasis constellation of H and A but 
admits any p(0), the other requires a generic eigenbasis 
constellation of p(0) and A but admits any H. 

In conclusion, von Neumann’s demonstration of ther¬ 
malization for isolated many-body systems has been gen¬ 
eralized to arbitrary observables. The remaining prereq¬ 
uisites for thermalization are thus rather weak, namely 
a Hamiltonian with generic eigenvalues and a generic 
orientation of its eigenvectors relative to those of A, 
while the initial state p{Q) may still be chosen arbitrarily 
(mixed or pure, far from equilibrium or not). The first 
requirement (regarding En) is by now well established 
[3-lll|. and further generalizations like in Refs. fT3 - [l5| 
seem possible. With the second requirement, von Neu¬ 
mann essentially anticipated the foundation of random 
matrix theory [3| , which is very difficult to justify rigor¬ 


ously, but is extremely successful in practice, and can be 
corroborated by various intuitively convincing arguments 
M- For instance: Since our mind is used to thinking 
about the physical world in terms of individual “par¬ 
ticles” , we mostly come up with single-particle observ¬ 
ables A or sums thereof (kinetic energy (temperature), 
density, pressure, magnetization, etc.), whose eigenvec¬ 
tors are thus single-particle product states. In contrast, 
a generic Hamiltonian H includes particle-particle inter¬ 
actions, giving rise to a “completely different” eigenbasis 
without any “special relation” to that of H [^. 

While von Neumann had in mind a preset H and a 
varying (or typical) eigenbasis of A [^, the mathemati¬ 
cally equivalent but physically opposite viewpoint (fixed 
A, varying H) was emphasized, e.g., in Refs. (3l.l4 ll6l . ll7l| . 
Here, both views have been merged and significantly gen¬ 
eralized by treating all three operators A, H, and p(0) 
on an equal footing: After selecting two of them and 
assuming they exhibit a typical eigenbasis constellation, 
we were able to draw conclusions which are then entirely 
independent of the third one. Along these lines, we estab¬ 
lished the general new results (ED and ED concerning 
the generic long-time behavior (equilibration) for expec¬ 
tation values of arbitrary (even untypical) observables 
and the entropy of the concomitant equilibrium states. 
Furthermore, our findings H14I) . (1151) significantly general¬ 
ize previously known typicality results for arbitrary (even 
time-dependent) Hamiltonians. As a by-product we thus 
obtained a unifying framework for several key aspects 
of thermalization, such as the validation of the eigen¬ 
state thermalization hy pot hesis [3, [13, [13, by means of 
random matrix theory 2l|, recent explorations of “equi¬ 


libration” iDU and “canonical typicality” [Isl [13, H 
[39i , the lon g-la sting misjudgment of von Neumann’s work 
[3J, Wa, HO, El , and its rehabilitation in Ref. [3| . 
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